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Foreword 



This report contains the notes for a series of lectures 
given in the Physics Department at Purdue University in the 
period tourch 25 to April 5, l-)57. The writer is very grateful 
for the hospitality which he received at Purdue during thesf 
lectures . 



Nuclear astrophysics is a very fast-developing field. 
These notes have been prepared primarily for the use of present 
and potential research workers in this field. They attempt to 
present a rev'ew of the field with reasonable technical com- 
pleteness and to indicate where further research is needed. 

//hile many new «esult3 are given here, this report Is in no 
sense intended to be a publication of them, particularly since 
many of the results are very preliminary and it is expected 
that they will be considerably improved in the near future. 
Nevertheless, several important conclusions about stellar 
evolution can oe drawn from the results in preliminary form, 
and it has been l'tlt useful to present conclusions at this time 
which will continue to be at least qualitatively correct. 

No attempt has been made to polish the presentation 
of material ’n this report. The material has been basically 
designed for presentation in the form of lectures. Hence 
there has been no effort to document fully every statement, 
although a reasonable guide to the literature may be found 
in the references which are quoted. 

The writer has benefitted greatly In his understanding 
of much of the material presented here from conversations and 
correspondence with many colleagues both at Chalk River and 
elsewhere. Particular thanks are due to J.L. Greenstein, 

M. Schwarz3chi Id, <V.A. Fowler, G.R. and E.M. Burbidge, 

L.G. Elliott, and .V.T. Sharp. 



A. G. W, Cameron 



Chalk River, Ontario, 
25 April, 1957. 
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1. INTRODUCTION 



In these lectures we will be Interested in the development 
of different kinds of stars, in the nuclear reactions which can go 
on in their interiors, and in the bearing of these considerations 
on the chemical composition of the universe and the origin of the 
elements, we shall see that a good case can be made for the theory 
that the elements have been and are being made in stellar interiors* 
However, we must first briefly survey certain fields of astronomical 
knowledge which will particularly concern us. 

1.1 Clusters of galaxies. 

The largest organized units of matter in spaco appear to 
be the clusters of galaxies. Most galaxies appear to belong to 
clusters. A cluster may contain only a few galaxies, but there is 
a continuous range of sizes extending up to clusters with a 
membership of many thousand galaxies. The galaxies in the larger 
clusters have a greater internal velocity dispersion than do 
those in small clusters. The clusters themselves do not have large 
relative motions. Zwicky (1956) has pointed out that there Is a 
limit to the size of organized material units given by 

D c l < 2. Dg , 

v g 

where D c ^ is the diameter of a cluster, Dg is the average 
intergalactic distance in a cluster, Vg is the mean velocity of 
a galaxy and c is the velocity of light. The largest clusters have 
diameters of the order of D<j4« 

1.2 galaxies . 



The most numerous class of galaxies is that of the 
elliptical galaxies. These appear to be composed mostly of very 
old stars, with only relatively small amounts of interstellar gas 
and dust also present. The ellipticals vary in shape from those 
which are nearly spherical to those which are quite flattened. The 
spiral galaxies are more massive than the ellipticals. They have, 
in addition to a subsystem of old stars resembling an elliptical 
galaxy, spiral arms closely associated with large quantities of 
gas and dust. The spiral arms contain also many newly formed stars. 
Some smaller irregular galaxies also exist; these contain mostly 
stars and interstellar material similar to that found in spiral 
arms. The mass of a typical galaxy is lOHsolar masses, but 
galactic masses range both up and down from this figure by large 
factors. The stars typically found in spiral arms are called 
Population I; those typically found in elliptical galaxies are 
called Population II stars. 
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1,3 Star clusters. 



There are major subsystems in galaxies called star 
clusters. These occur in two classes: the globular and galactic 

clusters. The globular clusters contain Population II stars. A 
typical cluster may contain 100,000 3 tars, but some of them 
contain many millions ol stars. These clusters are nearly 
spherically distributed about the center of our galaxy. The 
galactic clusters are concentrated toward the plane containing 
the spiral arms and they contain Population I stars. Globular 
clusters are so massive and so compact that they are very stable 
against disruption due to the gravitational effects of passing 
stars. On the other hand, galactic clusters can be dispersed 
in space by such encounters in the course of a few billion years, 
an extreme form of galactic cluster is the 0 or B association, 
which consists of a few usually bright stars receding from 
each other with quite high velocities. These stars were in the 
neighborhood of a common point in space a few million years 
ago, but in a few more millions of years they will be so dispersed 
that it will be difficult to recognize them as an associated group 
of s tars . 

1.4 Stellar Luminosities. 



There are many different classes of stars, and during 
the course of these lectures we shall try to reach an understanding 
of the paths of development of some of these stellar classes. 

One of the most remarkable properties of stars Is the wide range of 
their luminosities. The brightest stars known emit about 10® as 
much energy per second as the sun; the faintest known omit only 
10"® as much energy per second as the sun. The total rate of 
energy emission of a star Is called it3 absolute bolometric 
magnitude. Astronomical magnitudes form a logarithmic scale 
with a factor 2,512 between magnitude classes. Because of the 
limited spectral transmission of the earth’s atmosphere and the 
even more limited spectral sensitivity of most recording instruments, 
it Is difficult to translate measurements Into bolometric mag- 
nitudes. Hence in practice there are many systems of practical 
magnitudes which are based on specific spectral r anges , These 
magnitudes include the visual, photographic, and photoelectric 
magnitudes in various colors. The difference between the 
photographic and visual magnitudes of a star is called its color 
index. This gives a measure of the surface temperature of the 
star. The absolute magnitude of a star is equal by definition 
to the apparent magnitude it would have if it were placed 
10 parsecs away from us. The parsec Is a unit of astronomical 
distance equal to 3.258 light years. 
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1.5 Stellar spectra. 

The surface temperatures of the stars vary from slightly 
less than 2000°K to more than 500 # 000°K. The elements present in 
the surface layers of a star occur in various stages of ionization 
at different temperatures; hence the appearance of a stellar 
spectrum is sensitive to the surface temperature of the star. The 
spectra of most stars can be arranged in the following continuous 
sequence (each class is subdivided into tenths): 



Cla ss 
Class 



0 : 



Class A: 



Cla33 F: 



Class G: 



Class K: 



Temperatures of 25,000°K up. 
helium. 



Lines of ionized 



25,000 - 11,000°K. The lines of hyurogen and 
neutral helium are conspicuous at BO. Ionized 
oxygen and ionized carbon become strong at B3. 
Neutral helium is strongest at B5. Hydrogen 
lines become progressively stronger in the 
higher numbered subdivisions of the class. 

At AO hydrogen and ionized magnesium lines are 
strongest, while the helium and ionized oxygen 
lines have disappeared. Hydrogen lines weaken 
in the higher subdivisions, while ionized metals 
(Fe, Ti, Ca) strengthen. 10,700 - 7500°K. 

Class FO is rich 5n lines of the ionized metals, 
the strongest being the H and K lines of singly 
ionized Ca. Metallic lines strengthen and 
hydrogen lines weaken as we pass through this 
class. 7500 - 6000°K. 

In this class the lines of the neutral metals 
become strong while the hydrogen lines continue 
to weaken. Molecular bands of CN and CH appear. 
6000 - 4910°K. Our sun is class G2. 

In general molecular bands and lines of neutral 
metals become much stronger while the lines of 
hydrogen and ionized metals continue to weaken. 
At K5 the lines of TiO are weakly visible. 

4910 - 3500°K. 



Class M: The characteristic feature is the complex spectrum 

of molecular oxide bands,, of which the TiO bands 
are strongest. 5500 - 2200°K. 



Some stars do not fit into this sequence of spectral classes 
Therefore some additional spectral classes have been established which 
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parallel the above classes in temperature. These classes include 
the following: - 

Class S : A low temperature class parallel to Class M. 

This is still characterized by molecular oxide 
bands but the most prominent feature is the ZrO 
bands. Certain elements such as Zr, Y, Ba, La 
and Sr give strong atomic lines and oxide bands. 
Lines of neutral technetium are usually seen* 

Classes R and K (or Class C): Parallel in temperature 

to the ordinary classes K and M. The spectrum 
is characterized not by oxide but by molecular 
carbide bands, such a3 those of CN, Cg, and CH. 

Class W: Extremely high temperature objects called Volf- 

Rayet stars with bright, broad, hazy emission 
lines of ionized helium and highly ionized carbon, 
oxygen, and nitrogen. Two sequences exist: the 

WC stars have strong carbon lines and weak nitrogen 
lines; in the WN class the reverse is true. 

1.6 Surface temperature - luminosity diagrams. 

It has been known for a long time that if one plots some 
form of stellar magnitude against some measure of surface temperature, 
then most points tend to cluster along certain lines in the diagram. 
This type of diagram is often called a color-magnitude diagram or a 
Hertzsprung-Rus-sell (H-R) diagram. In the last few years it has been 
found that characteristically different H-R diagrams are obtained 
for globular and galactic clusters and for other types of associated 
objects. These will be discussed in greater detail below. However, 
at this point it may be stated that most stars cluster about a 
single line on these diagrams, called the main sequence, 

1.7 Mass -Luminosity relation. 

If two stars are joined to form a binary pair and if 
the periods, dimensions, and orientations of their orbits can be 
determined, then the masses of the stars follow from Kepler's 
third law. Quite a few stellar masses have been determined in this 
way for visual binary pairs. Russell and Moore (1940) have given 
us an empirical relationship between the mass and luminosity of 
main sequence stars. If the intrinsic luminosity L and the mass M 
are measured in solar units, then 

log L * 3.82 log M - 0,24 

The constant -0.24 expresses the fact that the sun is overluminous 
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for this relationship by 0.60 astronomical magnitudes. 

1.8 Stellar composition. 

There appears to be a striking parallelism between the 
compositions of the earth and the meteorites and the sun. This 
parallelism will be discussed in much greater detail below. 

However, at this point we may note that the earth arid meteorites 
have lost the volatile substances consisting of gases and substances 
with low boiling points such as mercury. This actually accounts 
for most of the mass, because it is found that the sun is nearly 
all hydrogen, with about 5 per cent of helium by numbers of atoms 
and much less than one per cent of other elements. Suppose we 
denote the atomic ratio of the metals to hydrogen in the sun by R 0 , 
One of the most important discoveries in modern astrophysics is 
that this ratio is not a constant of matter in cosmic proportions, 1 
Stars which can be classed as extreme Population II commonly have 
ratios of about R o /l0 or R o /20. In some extreme objects the ratio 
appears to be still much smaller. On the other hand, in stars which 
can be classed as extreme Population I the ratio is larger than in 
the sun, values of 2R 0 to 4R 0 being common. The light elements 
such as C, N, 0 and also heavier elements such as Ba appear to 
vary in abundance by about the same factors as the metals. 

1.9 Anomalous stellar abundances. 



Many stars have individual abundance anomalies which 
differ markedly from the solar abundances or the systematic 
Variations of solar abundances characteristic of extreme 
Populations I and II stars. The stars of spectral class S appear 
to have unusually large abundances of most elements heavier than 
about germanium. These heavy elements are also overabundant in 
a small group of stars called Ball stars (the designation "11“ 
means that the barium is singly ionized; "I 11 would have indicated 
neutral barium). Carbon is very overabundant in stars of spectral 
classes R and N, in "CH U stars, and to a lesser extent in Ball 
stars and in class 3. The "peculiar A" stars or spectrum variables 
have strong surface magnetic fields and an associated large 
overabundance of heavy elements and, in particular, of the rare 
earths. Wolf-Rayet stars (class M) appear to be deficient in 
hydrogen but rich in helium, carbon or nitrogen, oxygen, and 
often neon. Me shall be interested in seeing how these abundance 
anomalies can arise as a result of nuclear reactions which go 
on in stellar Interiors and at stellar surfaces. 

1.10 Interstellar matter. 



Something like half the mass of our galaxy exists in the 
form of gas and dust between the stars. Most of this gas is in the 
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form of atomic hydrogen. To the extent that the abundances 
of other elements can be determined, they arc of the same order of 
magnitude as in the sun. Tho interstellar matter is a very chemically 
reactive medium, and much of the heavier elements have formed 
chemical compounds which have collected to form dust grains. There 
is a considerable uncertainty as to the size of these grains. The 
interstellar absorption of star light is caused by the scattering 
properties of the dust grains. 

1.11 The ages of the sun and s tars . 

The presence of naturally radioactive material on the 
earth shows that the material from which it was formea does not have 
an infinite age but that there was a process in which the elements 
were 'formed. Considerations of the amount of energy which can be 
released in nuclear reactions allows limits to be placed on the 
length of time that stars can radiate energy at their presently 
observed rates. There are many dynamical properties of the stars in 
the galaxy which aro functions of a galactic time scale. These 
include the time required to disrupt galactic clusters ana binary 
stars. From these and many other considerations it appears that 
the Population II stars in our galaxy were formed about 7 x 10 9 years 
ago. The Population I stars have been forming continuously 3ince 
then. Our sun is an old Population I star. Our planetary system 
has an age of 4.5 x 10 y years, and another period of a t least 
0.5 x 10" years passed between the time of formation of the elements 
from which the solar system was formed and the time at which the 
formation process caused chemical separations of uranium and lead, 

Tho v sun probably formed at the 3ame time as the planetary system. 

The dynamical properties of 0 Associations show that their stars 
have been formed as recently as 10® years ago, 

1.12“ Ejection of mass from stars. 



We have seen that the stars are being continuously formed 
out of the interstellar medium. They aro also in the process of 
ejecting material to the interstellar medium. There are three 
principal mechanisms by which this can take place: 

Supernova explosions: At the peak of its light curve a 

supernova usually outshines the galaxy in which it is situated. 
Spectroscopic observations show that a substantial fraction of a 
solar mass is ejected in a supernova explosion. The ejected gases 
have outward velocities of some thousands of kilometers per second. 
About one supernova explosion per 500 years occurs in an average 
galaxy. 



Nova explosions: These are much leas spectacular than 

supernova explosions; only ahout 10"’ as much light is emitted, about 
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0.001 solar mass Is ejected in an outburst. Some of the novas which 
had smaller light amplitudes have been observed to flare up more than 
once, and it is likely that most novas are recurrent. The time 
between recurrences is of the order of a century, but there is a very 
large dispersion in this interval. 

Continuous ejection: after a star has evolved away from 

the main sequence there appear to be many stages In which clouds 
of gas can be emitted from the surface. Deutsch (1956) has observed 
a case in which a class M supergiant star is emitting matter at the 
rate of one solar mass per 30,000,000 years or less. Even white 
dwarf stars have been observed to emit matter (Greens tein 1956). 

This ejection appears to be associated with intense magnetic fields 
on the surfaces of the emitting objects. 

1.13 The formation of the elements. 



We have seen that stars are continually being born in 
interstellar 3pace and that in the later stages of their evolution 
they return matter to the interstellar medium. Galpoter (1955) has 
estimated that the amount of matter which has been formed into stars 
and has been returned to the interstellar medium is of the same 
order of magnitude as the total mass of all presently existing 
stars. We have also seen that the newly-formed 3 tars have much 
larger abundances of heavier elements than the old stars of 
Population II, and hence presumably the interstellar medium has 
been similarly enriched In heavier elements. During the course 
of these lectures we shall see that there are many types of 
nuclear reactions which can go on in stellar interiors which can 
produce certain of the elements with their observed relative 
abundances . 



It is a natural conjecture from the above considerations 
that the elements have been formed from hydrogen during the lif r 
of the galaxy. Thus about 7 x 10® years ago it is suggested that 
the galaxy was a large mass of gas (all or mostly hydrogen) in 
space. Much of this gas condensed to form Population II stars; 
the first stars to condense were very massive and evolved very 
quickly, ejecting heavier elements into the interstellar medium. 

These elements were Incorporated into other stars and further 
processed by nuclear reactions. The process of element formation < 

should therefore be still going on in the stars today. 

2. the development of the galaxies. 

Strictly speaking, the theories of element formation 
to be described are nearly independent of cosmological considerations 
because element formation is believed to take place after the formation 
of stars in the galaxy. This is compatible with the leading cos- 
mological theories in the present stage of development of the theories 
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of origin of the elements. However, there are some cosmology-linked 
effects that are susceptible to both observation and calculation, 
and it seems useful to review briefly some of the ideas about 
cosmology and the development of galaxies. ; 

2.1 Cosmology. 



To illustrate the main differences likely to be important 
to theories of element formation, we can briefly consider the 
cosmological exploding models (without cosmological constant) of 
Friedman (1922), and the steady state model of Bondi and Gold 
(1948) and Hoyle (1948, 1949). These models give different second- 
order terms in the expected red-shift-magnitude relation for 
distant galaxies. Present techniques of photoelectric photometry 
are promising to subject these second order differences to 
observational test. 

In the exploding models the universe was once in a 
highly compressed state. It was at one time thought (alpher and 
Herman 1950) that the conditions existing in the compressed state 
were suitable for formation of the elements In a non-equilibrium 
process of neutron capture. The difficulty with such theories is 
that there is no means of building nuclei heavier than He* unless 
one goes to unrealistically high densities of matter. Nevertheless, 
it may be that one should expect primordial abundances of deuterium. 
He 3 and He 4 to be formed In the early stages of an exploding universe. 
This can be tested In principle by measuring the helium abundances 
and deuterium abundances in Population II stars. The problem Is 
complicated because we know that helium Is one of the main products 
formed in stellar thermonuclear reactions, and helium abundances 
are, very difficult to determine accurately from spectroscopic 
measurements. Deuterium is even harder to measure ih the stars 
because its lines are shifted only slightly from those of hydrogen 
and its abundance is small compared to that of hydrogen. Perhaps 
we will have to wait for radio astronomy techniques to improve In 
sensitivity by another order of magnitude so that one can hope 
to detect the neutral deuterium in interstellar space. 

In the early stages of the exploding universe the density 
of radiation exceeds that of matter. Gamow (1953) has pointed 
out that no local condensations of matter are possible until the 
expansion has proceeded far enough for the radiation dens ity to 
fall below the matter density. Then local condensations are 
possible if the matter density is large enough. 

In the steady ;state model the average density of matter in 
space is constant in tinje, the expansion of the galaxies being 
balanced by the creatiori of new matter in space. In this picture the 
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new gas formed in intergalactic space, augmented by matter ejected 
from stars in small galaxies, condenses to form new galaxies. The 
gas in these new galaxies would have small abundances of all the 
elements present initially. These small abundances would have 
important effects in the stars forming in the new galaxy which 
may be tested by calculations in several ways. 

Preliminary photoelectric red-shift measurements of 
Baum (1957) show that the red-shift-magnitude relationship is 
very nearly linear out to velocities of recession of 0.4 c. 

These measurements support one of Friedman's models in which the 
expanding universe is closed; it is slowing in its expansion 
and will eventually stop and then start contracting. The data 
seem to be nearly good enough to disprove the steady state model. 

2.2 Fragmentation of gas clouds into galaxies and stars. 

We have seen that either of the cosmological theories 
is likely to produce a mass of gas, mostly hydrogen, which is 
unstable to a contraction into galaxies. Hoyle (1953) has given 
an attractive theory for the course of development of this gas. 

His theory is outlined briefly below. 



The key to Hoyle’s galactic evolutionary theory lies 
in the thermal properties of hydrogen gas. Let us consider a 
hydrogen gas cloud with an initial density f> 0 - 10~ 27 gm./cm. 3 
To estimate the temperature of the gas, we must realize that 
most cosmic objects have very large Reynold's numbers and hence 
the gas is likely to be turbulent, i.e. to develop mass motions. 

The energy of the mass motions becomes converted into heat. 

Hoyle estimates that certainly more than lO^- 2 ergs/gm. and possibly 
more than 10 14 ergs/gm. is made available in this way. 

The thermal energy of atomic hydrogen is £ RT per gram, 
where R is the gas constant and T is the temperature (°K). Hence 
somewhat more than 10^ 2 ergs/gm. is required to heat hydrogen from 
0°K to 10,000°K. At this point ionization of the hydrogen sets 
in. Ionization is nearly complete at 25,000°K. It requires some- 
what more than 10^3 ergs/gm. to ionize hydrogen. 



Hydrogen gas can radiate energy in three ways: (1), line 
emission, which will be neglected, (2), free-free transitions, in 
which unbound electrons make transitions between two states of the 
continuum, and (3), free-bound transitions, in which electrons are 
captured from a state of the continuum into a bound state of the 
hydrogen atom. The rate of hydrogen radiation from the latter two 
processes is 



1.45 x 1CT 27 T 1 / 2 \l+ 3.85 x 10 5 

T 






— 



N 



2 ergs/cmP sec.. 



\ 
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where N is the number of hydrogen atoms per cm. and T is the kinetic 
temperature in °K. This formula holds well only for T >25,000 °K. 

We see that the rate of radiation decreases until very large tempera- 
tures are reached. 

From these properties of hydrogen gas we can see that the 
large mass of hydrogen gas with which we have to deal is lively to 
have a temperature of either close to 10,000°K or greater than 10 5 °K. 
The reason for these two temperature classes is that the time of 
development of a galaxy is of the order of 10^' seconds, and in this 
time hydrogen Initially at temperatures -'10° °K can cool by 
radiation until most of the hydrogen atoms have recombined with 
electrons • 



Let us consider first the case of a gas cloud at a 
temperature of 10,000 °K. The necessary condition for a gas cloud 
of uniform composition to be in thermal equilibrium, in the case 
where radiation pressure is negligible so that gravitational forces 
are balanced by gradients in the gas pressure, is that the total 
gravitational potential energy JTX (taken as positive) mu3t be equal 
to twice the thermal energy of the gas. If is greater than this 
then a spherical cloud will contract. This gives us the following 
contraction condition: 



GM 



y 



3KT, 



( 2 . 2 . 1 ) 



where M i3 the total mass and V the total volume of the gas cloud and 
G is the universal gravitational constant. This inequality gives a 
lower limit to the mass of a gas of neutral hydrogen Q which can condense. 
For M = |0 o V * 10“27v, then (2.2,1) gives 3.6 x 10 y solar masses 
as the minimum mass of hydrogen gas which can condense. During 
such condensation the thermal energy is quickly radiated «way and 
the temperature remains at 10,000 °K. Hence the contraction of the 
gas takes place under isothermal conditions. 

The thermal energy released in a uniform isothermal 
contraction from a volume V to a volume V* is 



MRT 



A v 



( 2 . 2 . 2 ) 



The gravitation energy released by the contraction is of the 
order of 



GM 2 




1/3 



1 



1 



(2.2.3) 
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Now if the initial mass of the hydrogen gas cloud is large compared 
to the minimum of (2.2.1), then it follows that (2.2.3) is large 
compared to (2.2.2), regardless of the value of V*. Hence in the 
contraction of such a cloud only a small part of the gravitational 
potential energy is converted into thermal energy, and the rest 
must go into mass motions which will be adequate to expand the 
cloud back to nearly its original dimensions. Such a cloud 
therefore cannot undergo appreciable contraction. However, 3uch 
a cloud is also unstable to the fragmentation into subcondensations 
which have about the minimum mass consistent with (2.2.1). We 
may identify such a fragmentation process with the formation of 
a cluster of galaxies, the subcondensations having a mass 
--•3.6 x 10® solar masses as noted above. These are about the 
masses of elliptical galaxies. 

Let us now consider the condensation of a galaxy of 
this minimum mass. Expressions (2.2.2) and (2.2.3) remain 
comparable for V' = 0.1V; the former is still 30 percent of 
the latter for V' = 10"^V, and it has fallen to 5 percent of the 
latter for V’ = 10“®V. Hence there can be an almost complete 
dissipation of the energy of contraction for a decrease of linear 
dimensions of a factor 2 or 3, but further contraction is then 
prevented by the development of mass motions. However, after a 
contraction of linear dimensions by about a factor of 3, the 
galaxy is then unstable to a fragmentation into further 
subcondensations . These in turn can contract in linear dimensions 
by about a factor of 3, followed by further fragmentation, and 
30 on. as the density of the subcondensations increases, the 
rates of contraction and radiation also increase. If a sub- 
condensation contracts by a factor of k^/3, and then divides 
into k equal masses having a radius equal to l/k of that of 
the original subcondensation, and if many such stages of 
condensation take place, then Hoyle points out that the time for 
all stages of condensation Is, in units of that required for the 
first stage, 



1 + 3^ + 1. + _1 + ... = k 
k k 2 k 3 k-1 

For k 2 / 3 ^3, k^5, and the time required for an infinity of steps 
is only 25 per cent longer than that required for the first ste-p. 
Hoyle postulates that the condensation of the first step will take 
a time of the order of 10 17 seconds, and since the fragmentation 
will really be a continuous process, very small subunits will be 
formed over the whole of this time scale. 

The process of fragmentation will stop when the density 
of the subunits has become so large that they are opaque to their 
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own radiation. Hoyle shows that such final subunits have masses 
which generally should lie in the range 0.3 to 1.5 solar masses, 
with some dispersion which will ^extend outside this mass range. 
After a galaxy has been fragmented into a heirarchal structure 
of subcondensations, these must fall toward the galactic center 
of gravity where the heirarchal structure will be broken up by 
mutual interactions. This causes a shrinkage of the Dimensions 
of the galaxies. Occasionally certain subcondensations will 
resist disruption, and it is tempting to identify these with the 
globular clusters. 

Let us now consider a gas cloud which has a temperature 
initially in excess of 10° °K. The minimum mass required for 
condensation of this gas cloud is given by 

~ /_ > 6RT, (2.2.4) 

VV3 

where the altered form of the heat content refers to .fully ionized 
hydrogen. For M = 10“ 27 V the minimum mass is 5 x 10 11 solar masses* 
As contraction occurs the rate of radiation starts to exceed the 
rate of release of thermal energy, and eventually a point of 
temperature instability must be reached that swiftly reduces the 
temperature to 10,000 °K. Assuming that the linear dimensions have 
decreased by a factor of 5 before this occurs, then the first frag- 
mentation of this gas cloud will be into dwarf galaxies with masses 
of about 3 x 10® solar masses. Further fragmentation will then 
occur as described above. These dwarf galaxies will collide with 
one another and most will be disrupted; some will escape into 
intergalactic space. We may identify the massive spiral galaxies 
with a formation process of this type. 

Not all the gas in these galaxies will be formed into 
stars; there will be less tendency for condensation in regions of 
smaller gas density. Such gas will contract to the center of the 
galaxy, or, if the galaxy has a fairly large amount of total 
angular momentum, then the residual gas will contract to a disk 
passing through the center of the galaxy. Further star formation 
can then take place. 



3. HERTZS PRUNG-RUS S ELL DIAGRAMS . 



It has been remarked above that the stars are not uniformly 
spaced in a diagram in which their luminosity is plotted aga'nst their 
surface temperature. Many conclusions about the course of stellar 
evolution can be drawn from a study of these Hertzsprung-Russell 
diagrams for individual star clusters. The methods of plotting the 




the diagrams take a wide variety of forms : the luminosity can 

be expressed in bolometric, photographic, visual, or photoelectrio 
magnitudes in different colors . The surface temperature can be 
expressed as a temperature, spectral class, color index, or as a 
difference of photoelectric magnitudes in different colors. However, 
the basic star sequences in these different types of H-R diagrams are 
very similar, only relatively minor distortions being introduced by 
the different forms of presentation* 

3,1 H-R diagrams for galactic clusters. 

Figure 3,1,1 shows a composite diagram of several 
galactic clusters (Sandage 1956), The common line which these 
clusters share on the left is the main sequence. The short 
segments of these diagrams on the upper right lie in the red 
giant and red supergiant region. It may be noted that all the 
galactic clusters except M67 have a gap between the sequences 
of stars lying on or near the main sequence and their red giant 
branches. This is called the Hertzsprung gap. The gap is wedge- 
shaped, being non-existent for M67 and becoming larger as one 
goes upwards in the diagram. 

Many features of stellar evolution may be deduced 
from an inspection of Fig. 3.1.1 The stars at the top of the 
diagram are using up their energy sources most rapidly. Hence 
the clusters at the top of the diagram must have been formed 
very recently. We can see that the brighter stars in these - 
clusters lie on sequences which curve away from the main 
sequence upwards and to the right. This establishes the 
direction of evolution of the individual stars when they evolve 
away from the main sequence. After such stars reach the boundary 
of the Hertzsprung gap they must evolve to the right very quickly 
to reach the giant branch appropriate to the cluster. They then 
sit on the giant branch for a considerable period of time. It 
may be noted that In most of the clusters the giant stars have 
about the same luminosity as those which are evolving off the 
main sequence. However, in an old galactic cluster like M67 
the stars brighten by about three magnitudes after evolving away 
from the main sequence. The lengths of time one can expect 
the stars to remain on the main sequence are shown at the right 
of Fig. 3,1.1, The stars in galactic clusters belong to Population I 
those In M67 can be described as old Population I, while those in 
clusters such as NGC2362 and h +-X Persei have been formed very 
recently and are extreme Population I. 

Our composite diagram does not suffice to tell us 
anything about the tracks of the 3 -cars before they reach the 
main sequence or after they have exhausted their sources of 
energy on the giant branch. Since we expect that the stars form 
from gas clouds contracting from Interstellar space, then initially 
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Figure 3.1.1 




a composite color-magnitude diagram for galactic clusters* 
In this diagram visual magnitudes are plotted against the 
difference of the photoelectric magnitudes measured In 
blue and visual light. For comparison, the giant branch 
of the globular cluster M3 Is also included. The lengths 
of time required for stars to evolve away from the main 
sequence are shown on the right. 
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such contracting gas clouds must have very low surface temperatures 
owing to the large surface area available to radiate the energy 
released in the contraction. This would place the forming objects 
far to the right of the main sequence in the H-R diagram, as 'the 
contraction proceeds the stars move towards the left until they 
reach the main sequence. Since the more massive stars radiate 
energy enormously faster than the less massive ones, we can expect 
them to reach the main sequence first. This raises the possibility 
that one might be able to find galactic clusters so recently formed 
that their brightest members would lie on the main sequence and the 
faintest members would lie to the right of the main sequence, 
talker (1356, 1957) has found several such clusters. The color- 
magnitude diagram found by him for NGC2264 is shown in Pig. <5.1.2. 
The brightest stars lie on the main sequence; it may be noted that 
the five giant stars are here less luminous than the brightest main 
sequence stars. The fainter stars lie above and to the right of 
the main sequence. Many of these stars show light variability and 
bright hydrogen emission lines. This places them in the class of 
"T-Tauri variables", which are believed to be stars in gravitational 
contraction which are still interacting with the interstellar medium 
from which they were formed. 

3,2 H-R diagrams for globular clusters. 

There are many characteristic individual differences in 
the color-magnitude diagrams for the globular clusters. These 
clusters contain Population II stars formed a long time ago; like 
the galactic cluster M67 they do not show a Hertzsprung gap. M67 
is compared with the globular clusters Mu and M92 in Pig. 3,2.1. 

The main sequences of the globular clusters are not shown, but 
it is known that they lie to the left and below the main sequence 
of the galactic clusters (it appears that the main sequences of 
different globular clusters lie in somewhat different positions in 
the H-R diagram). Stars evolve away from the main sequence at 
about the same luminosity as in M67, but they increase in luminosity 
to a much greater extent when they reach the red giant region. It 
can be seen from Pig. 3.2.1 that after reaching the tip of the 
giant branch at the upper right, the stars reverse their direction 
of evolution, decrease somewhat In luminosity, and move to left 

in the diagram. This region of the diagram is called the horizontal 
branch. The horizontal branch contains a region, shown by the gap 
in the M5 diagram, where the stars are light variables; such stars 
are called cluster variables, 

Baum(1957) points out that the features of the color- 
magnitude diagram that differ from one globular cluster to another 
include (1), the position of the main sequence, (2), the color of 
the top of the red giant branch, (3) the shape and tilt of the red 
giant branch, (4), the number of cluster type variables, (5), the 
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Figure 3,1.2 




b-v 



Color-magnitude diagram of NGC2264. Here photoelectric visual 
magnitudes are plotted against the difference of blue and 
visual magnitudes. The stars are reddened by interstellar 
absorption. Dots represent photoelectric and circles photo- 
graphic observations. Vertical bars show known light-variables 
and horizontal bars show stars with bright emission lines. The 
lines snow standard main sequence and giant branches for 
Population I stars. 
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COLOR IRDEX n-V 



Color-magnitude diagrams of the galactic cluster M67 compared 
with the globular clusters M3 and M92, The absolute visual 
magnitudes are plotted against the difference of the photo- 
electric blue and visual magnitudes. 
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populousneaa of the horizontal branch at the right of the cluster- 
type variables, and (6), the ahape and tilt of the horizontal branch. 
These differencea appear to be connected with differences in chemical 
composition between the atara in different globular clusters, and 
there are corresponding correlations in the differences in these 
features in different clusters. The abundance .ratio of metals to 
hydrogen varies quite a bit between different globular clusters but 
is always much less than in the sun, 

* 

Greenstein (1956) has made extensive jstudies of the spectra 
of subdwarf and white dwarf stars. The subdwarf stars are members 
of extreme Population II which lie below the Population 1 main 
sequence. The H-R diagram for these stars is shown in Pig. 3,2,2, 

It may be seen that the stars with the smallest abundance ratios 
of metals to hydrogen generally lie farthest below the Population I 
main sequence, 

3.3 Evolutionary tracks in the H-R diagram. 

Prom an inspection of the Hertzsprung-Russell diagrams for 
gldbular and galactic clusters we have reached certain conclusions 
about the evolutionary tracks of the stars in these diagrams, 

Sandage (1954) has obtained evolutionary tracks for globular cluster 
stars by mapping stars from their original positions on the main 
sequence to their presently observed positions. These evolutionary 
tracks are shown in Pig. 3,3,1, They confirm the” prevl ous discussion. 
It should be noted that the stars in the giant branch have brightened 
so much that their evolution is very rapid. Such stars have nearly 
the same mass. 

After stars travel to the left on the horizontal branch 
they are often quite unstable, shed mass, and sjhow peculiar 
abundances of certain elements, Pre- and post*novas exist in the 
far left region of the H-R diagram; the nova explosions appear to 
be connected with Internal instabilities whloh^arlse near the 
end of their course of evolution. It appears ihat the stars finally 
take the form of white dwarfs at the end of th^. r course of evolution. 
Such stars are composed of a degenerate electron gas. They have a 
very high density and contain no hydrogen or energy sources in the 
interior. Their internal heat energy is gradually radiated into 
space and they are expected to cool off at constant radius. Luyten 
(1952) has plotted a color-magnitude diagram of the white dwarf 
stars; this is shown in Pig, 3,3,2, The white dwarf stars have a 
very low luminosity, but their small diameters imply a small surface 
area and hence they are usually blue objects. It may be seen that 
most white dwarfs appear to have diameters varying in size between 
one and four times that of the earth. The white dwarf stars are 
very numerous in space, second only to the numbers of stars on the 
main sequences. They represent the graveyard of stars which have 
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Figure 3.2.2 




Hertzsprung-Russell diagram of the subdwarf stars. The 
photographic magnitude Is plotted against the spectral 
class. The solid dots show stars with a smaller abundance 
ratio of metals to hydrogen than the open triangles. The 
line shows the position of the main sequence for Population 
I 8 tars. 
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Figure 5,5,1 




Evolutionary tracks In the Hertzsprung-Russell diagram for 
globular clusters. Absolute bolometric magnitudes are 
plotted against stellar surface temperatures » The heavy 
solid lines show the positions of the original main sequence 
and also of the present locus of stars in the globular 
cluster M3, 



21 



CRL-41 



Figure 3.3.2 




Col or -magnitude diagram of the white dwarf stars. The 
photographic magnitudes are plotted against the color 
Index. The thick line on the upper right sho*s the 
Population I main sequence. The other lines give the 
diameter of the white dwarf star in units of that of 
the earth. The track of a star cooling at constant 
radius would lie parallel to these lines. 
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evolved off the main sequence and exhausted all their energy sources. 
4. PHYSICAL CONDITIONS IN STELLAR INTERIORS, 



Efforts have been made for many years to determine the 
physical conditions at various points in stellar interiors. These 
attempts have improved with the passage of t J me and tne more realistic 
inclusion of various physical effects in the theory. The development 
of the subject has been governed by the fact that the differential 
equations describing the stellar interior are non-linear. Thus sane 
elegant mathematical analyses have been developed which enable stellar 
models to be computed using analytic approximations to certain physical 
properties; these models have the nice property that computations for 
one assumption of mass and luminosity can be transformed to other 
values of mass and luminosity by using certain dimensionless variables 
in the calculations. These methods have given very useful results, 
but they will not be described here because stellar models are now 
being calculated numerically on electronic computers, and different 
methods of analysis are best suited to computer calculations. 

4.1 Differential equations of the stellar interior. 

In most main sequence stars atoms in the interior are 
largely stripped of electrons and hence the gas obeys perfect gas 
laws except under conditions of very great density (> 10 4 gm./cm. 5 ) 
where electron degeneracy sets in. In order to write down the 
differential equations for the stellar interior, let us consider a 
thin shell of radius r about the center of the star with thickness dr. 
Then we get a condition for hydrostatic equilibrium In this layer 
by stating that the difference of pressure between top and bottom 
of the layer must support the weight of the layer: 




total pressure 

gas pressure p + radiation pressure p r , 
density in layer 
mass inside layer, 

universal constant of gravitation. 

The conservation of mass In the star is given by the 

relation: 



where P = 



P = 
M(r) = 
and G = 



dM(r ) = 4 ttt 2 p 

dr 



! 
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The conservation of energy in the star is given by the 

relation: 



dr 

where L(r) is net outward flux of energy entering the bottom of 
the layer and E is the energy production per gram per second in 
the layer. 



We will obtain another differential equation which 
describes the transport of energy in the stellar interior. In 
general there are three methods by which energy transport can 
take place: conduction, radiation, or convection. Conduction 

is not important unless electron degeneracy is present, a case 
that will be considered later. Let us consider the case where 
energy transport is solely by radiation. 

Energy transport by radiation takes place more rapidly 
for large temperature gradients than for small. It is impeded by 
high densities and high opacities of the material through which it 
must pass. These conditions are expressed in an equation for the 
temperature gradient required to transport a stated flux of energy 
(Chandrasekhar 1939) : 




where T is the temperature (°K), 

a is Stefan's radiation constant, 
e is the velocity of light, 

and K is the opacity of the material, which is an appropriate 
mean of the mass absorption coefficients of the materials present. 
This is discussed in more detail below. 



Now let us consider when energy transport by convection 
will take place. Let us consider a volume element of the gas which 
we will take large enough so that adiabatic conditions will hold if 
we give this element a sudden upward displacement. In this dis- 
placement pressure equi librium will be maintained. Under adiabatic 
conditions 



o 



P 1 ^ 



» 



where y = = ratio of specific heats. 

If the new density is less than that of the surrounding matter, the 
element will continue to rise, and hence that region of the star is 
unstable to the onset of convection. The condition for stability 
against convection is that the actual pressure gradient in the star 
must be less than that corresponding to a density gradient in 
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adiabatic equilibrium. In general, such stability requires 
relatively small temperature gradients. 

Under adiabatic conditions the temperature varies with 
the pressure according to the relation 

p(r-i)/Y. 



A rising volume element has a smaller density than its surroundings 
and hence it has a temperature excess. At higher levels it loses 
its excess heat by viscous effects and by radiation. Similarly, 
a descending element has too low a temperature and is heated in 
lower levels. This energy transport process will very quickly 
reduce a very steep temperature gradient until it is close to the 
limiting adiabatic gradient. Hence it may be seen that for steep 
temperature gradients, nearly all transport of energy is by con- 
vection; while for temperature gradients lower than the adiabatic 
value the energy transport is entirely by radiation. 

Thus for a case of convective equilibrium we have the 
adiabatic condition holding in the star to a good degree of 
approximation: 



1 dP = JC dp 
f dr o dr 

In general we can expect convection to be important 
near the centers of stars where thermonuclear energy generation, 
which often varies as a high power of the temperature, tries to 
establish a steep temperature gradient. Farther away from the 
center of a star where energy generation is not important the 
temperature gradient becomes less and energy transport is by 
radiation. However, convection can also set in in a region where 
a major constituent of the gas is partly ionized, as a volume 
element of the gas rises in such a region recombination takes 
place. This heats the volume element which continues to rise 
until the major constituent has become nearly all recombined. 

This process is responsible for establishing outer convection 
zones in stars where hydrogen is partly ionized. 

4.2 Additional equations for the 3tellar interior. 

In order to determine the structure of a star we need 
also to know some additional properties of the matter contained 
in the interior. One property is the equation of state of a 
perfect gas : 





T, 
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where k = Boltzmann's constant, 

H = mass of the hydrogen atom, 

and |i = mean molecular weight of- the matter, which is the inverse 
of the number of particles per unit atomic mass (p, = l/2 for hydrogen, 
4/3 for helium, and 2 for heavy elements). 



The radiation pressure is given by 



P r = 1 a T . 

3 

The stellar opacity K is a very complicated function. In 
an accurate calculation it is necessary to take into account as 
detailed functions of photon frequency all the scattering and 
absorbing processes which can take place in all the different kinds 
of atoms present. We shall not give a detailed theory of opacity 
but will content ourselves with mentioning the important physical 
processes which contribute to it. 



Under the conditions usually present in the interiors of 
main sequence stars the most important absorption process is 
photoelectric absorption in the heavier atoms present. Hydrogen 
and helium do not make appreciable contributions to the opacity 
in this process, since they are completely ionized. 

A lesser source of opacity is the free-free transition. 
In which an electron moving In the continuum Is excited to a 
higher state of the continuum. 



At higher temperatures an important source of opacity 
becomes photon scattering by free electrons (Thomson scattering or 
the Compton effect). 



Much work has been done with an approximation to the 
stellar opacity known as Kramers' formula: 

K = constant 

Here the constant is proportional to the abundances of elements 
heavier than helium. 

It may therefore be seen that the presence of heavier 
atoms In stellar interiors is very important from the point of 
view of the opacity. Higher abundances of heavy elements increase 
the opacity, which in turn raises the temperature gradients in the 
Interior. The rate of energy generation at the center is therefore 
Increased, and a star of given mass will be more luminous than 
one with a smaller content of heavy elements. We see here a reason 
why the main sequences of the globular clusters lie below that for 
Population I stars. 



I 
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In addition to a good opacity relation it is also 
necessary to obtain from physical theory the dependence of the 
energy generation rate E on temperature, density, and chemical 
composition of the material, as well a3 any contributions 
due to gravitational contraction. 

The four differential equations which apply at any 
point in a stellar interior, taken together with the additional 
relationships discussed here, allow a complete stellar model to 
be constructed, given only the total mass and distribution of 
chemical composition in the star. In practice such computations 
are very complex to ca^ry out, and really good results in which 
stellar evolution has been followed through a series of models 
with changing chemical composition are only now starting to be 
obtained with electronic computers. 

4,3 Electron degeneracy. 

The preceding discussion was based on the assumption 
that the perfect gas laws held to a good approximation in stellar 
interiors. This assumption fails at high densities owing to the 
onset of electron degeneracy. 

Let us consider electrons which are confined to a unit 
volume. According to quantum theory there are a limited number of 
states available to these electrons. The number of states 
available is 




in the momentum range p to p + dp. This allows two electrons to 
exist in each cell of phase space, one with each of the two possible 
spin orientations. The actual number of electrons in this momentum 
range will be less than or equal to the above quantity: 

N(p)dp ^ 12 ') 47C p 2 dp. 

W>) 

Let us compare this with the number of electrons in the same 
momentum interval which we would have for a Maxwell distribution 
of velocities: 



H( P )d P ‘ !• e -p 2 /emkT ^- p 2 dp> 

(2 7tmkT) 3 / 2 



where N e is the number of electrons per unit volume and m is the 
mass of an electron. Now at low momenta the Maxwell distribution 



- 27 



CRL-41 

violates the Pauli principle if 

j . y (4.3.1) 

(2rmkT) 3 / 2 h 3 

If N 0 = P , where p, e is the mean molecular weight per 

TeB- 

electron, then (4.3.1) can be expressed in the form: 

.o > 2(277'mkT) 5 / 2 = 8#1 x 10 -9 ^ e T 3 / 2 (4.3.2) 

h 3 



For p, e = 2 and T = 10® °K, the Pauli principle would be 
violated for densities in excess of 1.62 x 10 4 gm./cm.® 

Serious violations would certainly occur at the much higher 
densities sometimes encountered In stellar interiors. 

Therefore at high densities the Maxwell distribution 
is depleted of electrons with low momenta and enriched with 
electrons of higher momenta. At extremely high densities a 
good approximation to complete degeneracy must set in, in which: 



n(p) ws \ 4-rrp 2 , 

and N(p) = 0, 



(P< Po) 
(p>Po) 



(4.3.3) 



in which p Q is the Fermi threshold of the distribution 
(sometimes called the Fermi level, or the corresponding energy 
p 0 2 /2m is called the Fermi energy) . It is given by 



Po 



= r 3h£ 

l 8 < 




1/3 



Under conditions of interest in stellar interiors 
electron degeneracy sets In at high densities under conditions 
such that the electrons move with velocities which are a large 
fraction of the velocity of light. Hence relativistic mechanics 
must be used in a derivation of the equation of state of a 
degenerate electron gas. This has been done by Chandrasekhar 
(1939, 1901). He obtains the following results for the pressure 
and density of a degenerate .electron gas: 
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p = 7 Tm 4 c 5 f (x), (4.5.4) 

5h 3 

where f(x) = x(2x 2 - ‘6)(x 2 +D 1 / 2 + 3 sinh"^- x, 

and x = p 0 ; 

me 

= QlL » 5 e 5 - M-eHx 3 , (4.3.5) 

' 3h 3 



(4.3.4) and (4.3.5) represent, parametrically, the equation of 
state of a highly degenerate electron ga3 . The pressure In this 
gas is much larger than in a gas obeying the perfect gas laws. 

Such aadltional pressure will be referred to for convenience as 
electron degeneracy pressure. 

Highly degenerate matter in equilibrium under its own 
gravitat'on has some very interesting properties. The temperature 
of the matter is irrelevant as long as it is not so high that it 
destroys the good approximation to complete degeneracy. However, 
the degenerate matter is a good thermal conductor, so that the 
temperature distribution in the matter must be very flat. It turns 
out that there is an upper limit to the mass of such a completely 
degenerate body; the degeneracy pressure can support the weight of 
the ions and electrons provided that weight is not too large. The 
radius of the degenerate body approaches zero as this limiting mass 
is approached. The limiting mass is: 

5.756 solar masses. 



oh, for p e =2, M^1.44 solar masses. There are some corrections 
which should be made which somewhat reduce the value of the 
limiting mass given above. White dwarf stars appear to be highly 
degenerate bodies of the sort described to a good degree of 
approximation. The white dwarfs contain no energy sources, and 
gradually cool off at constant radius as their internal heat content 
is radiated to space from the thin non-degenerate atmosphere. 

One further property of degenerate matter should be 
mentioned owing to its importance in the subsequent discussion. 

Let us consider a volume element of the degenerate matter and let 
this volume element be suddenly compressed. The Fermi threshold 
of the electrons is raised at the higher density which results 



/ 
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from the compression. The internal energy stored in the electrons Is 
increased, and this energy can only come at the expense of the kinetic 
energy of the ions present which are not degenerate. Accordingly, 
the ionic temperature is reduced during a compression. Similarly, 
when a degenerate gas expands the ionic temperature increases. This 
behavior is opposite to that of an ordinary gas. If the temperature 
of the ions is high enough for the ions to take part in thermo- 
nuclear reactions, then the degenerate gas is very unstable. 



s. thermonuclear reaction rates. 



Although we have seen that electrons can form a degenerate 
gas In stellar interiors, we do not expect that the ions will become 
degenerate, owing to their much larger masses. Therefore, the ions 
will have a Maxwellian distribution of velocities under all conditions 
of interest. The tail of the Maxwell distribution decreases in an 
exponential manner, but nevertheless when considering thermonuclear 
reactions we are particularly interested in velocities in the far 
tail of the distribution, corresponding to particle energies often 
5 to 10 times kT. The reason for this is that nuclear reactions take 
place only when the colliding charged particles can penetrate their 
mutual Coulomb barrier far enough for the particles to come within 
the range of nuclear forces of each other. The probability for this 
penetration varies exponentially as the inverse of the velocity. 

The product of these two exponential factors has a sharp peak in the 
far tail of the Maxwell distribution. 

5.1 Nonresonant thermonuclear reactions. 



Most nuclear reactions have res onances; at certain Values of 
the bombarding energy the reaction cross section is much higher than 
at most other bombarding energies. We shall first be concerned 
with the case in which the thermonuclear bombardment energy region 
is far from the nearest resonance. However, the thermonuclear 
reaction rates will be formulated generally (Thompson 1957). 

The reaction rate r Is equal to the product of the number 
of possible pairs of reacting particles and to the probability that 
a reaction will take place between any reacting pair. For reactions 
between unlike particles it Is 

r = n x n 2 

and for reactions between like 

r = 1 n 2 

2 



particles It is 
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3 

Here n^, no, and n are the numbers of reacting particles per cm . 
The quantity<^d £> is an appropriate average of the product of the 
reaction cross section and the relative velocity of the reacting 
particles. It depends only on the temperature of the gas and the 
nuclear properties of the reaction. 



Let the particles have velocity distribution functions 
nj^v) and n 2 (v). Then 

r ■ jv? - vj | a (vj - v^) d®v^d^Vg # 

where the volume integrals cover both velocity spaces. Putting in 
Maxwellian velocity distributions, we obtain 



r = 



n l n 2 



m l 



,3/2 




d 3 ^ d 3 ^ 



Vi - v r 



<5 (vj - r 2 ) x 



2 2 

x exp - /l/2 m^v-^ + 1/2 m2V 2 



(- 



kT 



) 



where m^ and m^ are the masses of the two particles. 



We reduce the integral by changing to the relative velocity 
^ 3 v^ - ^2 and to the center of mass velocity 

V “ “i 7 ! * ’VZ 
* ®2 

The Integral over d V and over the solid angle in d v can be carried 
out Immediately. This gives 



r 



n l°2 4 IT / 



m l m 2 






3/2 



2-ff'kT (mi + m 2 ) J 



J. 



dv v 3 



a (v) exp 




^ . (5.1.1) 
kT 



At low bombarding energies and far away from resonances 
nuclear reaction cross sections can be expressed In the forms 



a * S exp - 2lf e^Z^Zg , 
E tiv 



where E is the relative energy which corresponds 
e is the electronic charge and Z^ and Zg are the 



to the 
atomic 



velocity v 
numbers of 



9 




t 
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the cojl iiuing particles. We write E = l/2 mv 2 , where m = m -, m 

m l +m 2 

S Is a slowly varying function of the bombarding energy which is 
considered to be constant in the present formulation of the problem. 
The reaction rate becomes : 

r = 2n-.no/ 2 \ l/g 3 fdE exp -/ E + B \ , 

llfm) (kT ) Jo ikT 

O 

where B _ 2/fZjZf>e , The integrand is sharply peeked and may 

i ^ 2 Tv 

be approximated by a Gaussian function. When this is done we get 
the following good approximation: 

fdx exp - (x + ax -1 / 2 ) ~ 2( 1 if )^ 2 ( 1 a) 1 ^ 3 exp - 3( 1 a) 2 ^ 3 

3 2 2 

This completes the derivation of the thermonuclear reaction 
rate except for some rearrangement to make the formulas convenient to 
use. We express the number density of particles in terms of the gas 
density by the relation 



n l = N oP x l # 
A 1 



where N 0 is Avogadro's number, A^ is the mass number of the particle, 
and x-^ is the concentration by weight of the particles of type 1 in 
the gas. We shall also express our energies in the center of mass 
system of the bombarding particles. Then p, the number of reactions 
per second per nucleus of type 2, is 



p = 434 p xi (Ax * A 2 > S T' 

A 



Ai 2 *2 Zj. z 2 



e sec . 



(5.1.2) 



where "T = 42.48 



Z] 2 Z 2 2 A x Ag"^ 3 



Ai + A2 




and 



S s SE exp 



988 Z Z 

X Cr 



-l. tz 

(Ai + A2)E 




(5.1.3) 



Here <f is the cross section for the reaction in barns. E is the 
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bombarding energy in the center of mass system in electron volts, and 
the temperature T is in units of 10° °K. 

The quantity S can sometimes be determined for a specific 
reaction by measuring the reaction cross section at the lowest possible 
bombarding energy and putting the measured value in equation (5.1.3). 

In other cases S must be calculated from resonance theory or from 
the statistical properties of nuclear resonances. This is discussed 
in greater detail below. 



Most reactions take place in the vicinity of the energy 
about which the approximated integrand is sharply peaked. This 
peak occurs at an energy 

= 1220 



2 2 c 

Zi Z 2 A 1 A 2 T- 



1/3 



electron volts. 



>»1 + Ag 



The full width at half maximum of the integrand is 



AE = 15.2 (TE^) 1 / 2 electron volts. 

This peak In the integrand will be called the "Gamow peak”. 



5.2 Resonant thermonuclear reactions . 

We now consider the case where there is a nuclear resonance 
close to the Gamow peak. In most cases of this sort in which we shall 
be interested the width of the resonance is very small compared to the 
width of the Gamow peak (here and in general the term "width” will be 
taken to mean the full width at half-maximum of a peek ) . In this 
case equation (5.1.1) can be written: 



r 



ait 



n i”2 / m i"fe 

1 2 rr kT l m x -Hflg ) 



3/2 




where v r is the relative velocity of the particles corresponding to 
the resonance. We see that the integral now is just the integrated 
cross section under the nuclear resonance. 



In a very refined treatment nuclear cross sections would be 
expressed in terms of the Breit-Wigner many-level formula. However, 
for most as trophysical purposes It is good enough to use the Breit- 
Wigner single level formula, particularly in the vicinity of a resonance 
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peak such as we are interested in here. The Breit-Wigner formula 
for the cross section of a reaction in which a particle a enters the 
nucleus and a particle b emerges from it is : 

a (a,b) = . S , (5.2.1) 

(E r -E > 2 + r 2 /i 

where \ a is the de Broglie wave length of particle a divided by 
2*7T , To# Tb» and Tare the partial widths of particles a and b and 
the total width of the level, r espectively, E is the bombarding 
energy, and E is the energy of the resonance, g is a statistical 
factor given r by 

g = 2J + 1 , 

(2s+l) (21+1) 



where s is the 3pin of the bombarding particle, I is the 3pin of the 
target nucleus, and J is the spin of the resonant level in the com- 
pound nucleus. 



by 



The integrated cross section is given to a good approximation 



r* (a.b) dE a ? rt z ^.f g Ta Tb , 

r 

where A must be evaluated at E^,, 



rfe write a subscript 1 for particle a and a subscript y 
for particle b (to indicate the fact that the resonance width is not 
likely to be smaller than the width of the Gamow peak unless the 
emergent particle is a photon). Then the resonant thermonuclear 
reaction rate becomes, where again the temperature T is in units 
of 10 6 °Kf 



p * 4.81 x 10 9 px 



A lT 3/2 



JjJjE Al + *2 \ 3//2 exp ( -0.0116 E r \ sec." 1 

Y Ul ±2 



(5.2.2) 



is the resonance bombarding energy in electron volts in the 
nter of mass system. 
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Resonant thermonuclear reaction rates are usually higher 
than nonresonant rates by a factor of something like 10°. This 
indicates the extreme importance of knowing whether there are 
resonances in the thermonuclear energy region when considering the 
as trophysical importance of specific nuclear reactions. 



r-r. 



ih? 



P. 

others 



In most reactions of interest the total width 
Usually one of these two partial widths is much smaller than 
In such a case the ratio fy/l ' can be replaced by the smaller partial 
width, and the larger partial width need not be known at all. 



5.3 Application of the Wlgner-Els enbud theory. 



Unless nuclear cross sections can be measured in the laboratory 
at energies close to the Gamow peak, and unless the closest resonance 
level is far from the Gamow peak, then the cross section at the Gamow 
peak must be determined using nuclear theory. A refined formulation 
of nuclear reaction theory has been given by V/igner and Sisenbud (1947), 
often called the dispersion theory of nuclear reactions. According 
to dispersion theory the partial width of a particle can be written 
in the form 



r 



= 2kPy 2 , 



where k is the wave number of the incident particle. 



k = 2.19 x 10 



9 



f A ^ Ag E 
+ A 2 



cm.“^ 



(5.3.1) 



(5.3.2) 



P is the barrier penetrability of the particle, and I® the reduced 
width of the particle in the level under consideration. In this theory 
the configuration space is divided into an internal region, to which 
all specifically nuclear effects are restricted, and an external region 
in which the system is composed of two particles interacting only 
through their Coulomb fields.' The barrier penetrability is 



P = 1 , (5.3.3) 

p2 + G 2 

where P and G are the regular and irregular solutions of the 
Schroedinger equation for a particle in a Coulomb field. For the 
low energies which are involved here, 

G 2 >> F 2 . 

The penetrability can be evaluated using the low energy 
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approximation given by Yost, Wheeler and Breit (1936), It is 
necessary to use two functions of the bombarding energy, measured 
as before in electron volts r 

p = kR, (5,3.4) 

and -\ = 157.4 Z-jZg AiA 9 X 1 ^ 2 (5.3,5) 

^ (A]_+Ag)E I 

Here R is the radius of the boundary between the internal and 
external regions in configuration space. In principle it should 
not matter what value is used for the radius of this boundary 
surface, provided it is taken large enough to avoid nuclear effects 
in the external region. However, in practice it is desirable to 
use as small a value as possible for the nuclear interaction 
distance, and it is customary to take it equal to the sum of the 
radii of the two interacting particles s 

R a 1.45 x 10" 13 (A^ 3 + Ag 1 / 3 ) cm. 



As we have seen, the p barrler penetrability for 
with angular momentum L is g£". This can oe written in 



part] cles 
the form 





'Z 



» 



where 1 


= 2 2 ^ 


(L 2 + A 2 )([L-1] 2 *\) ... (I 2 +'l 2 ) 




2 




i pmp 


V w \ l 

_ 




w- — 



We write q = 2^ j~(L^ + -n )^ ( [L-l] + \ ) »..(l^+ r ^)j« 

[(2L) i] S 



Uhder the conditions of interest exp 2rT-n s is large compared to unity 
and hence we can write to a good approximation 

2L 

exp - 





Por L * 0, q a l, We may evaluate the function 0^ in the low 
energy limit which is 



9 



L 





2L+1 



^2L+1 
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where x = (8 p^\ )^/ 2 an( j K 2 l+i(x) is the modified Bessel function of 
the second kind, 

Teichman and Wigner (1952) have found a practical sum rule 
limit for the reduced width y 2 whose order of magnitude is (dropping 
a factor 1.5 from Teichman and Wigner’s formula) : 




m R 



(5.3.6) 



where as before m is the reduced mass of the colliding particles 
and R is the radius of the boundary between the internal and external 
regions of configuration space. Lane (1954) has calculated reduced 
widths for nucleons in single particle states in light nuclei and 
finds values close to the upper limit given by the above inequality. 
Therefore, h”/mR will be called the single particle limit of the 
reduced widths. 



According to a modern view of nuclear level structure, a 
single nucleon with a given angular momentum can interact with a 
nucleus to form a single particle state about once in every 15 Mev. 
o ( f excitation energy. In practice the wave functions for these 
single particle states mix into the wave functions of all the nuclear 
states in the vicinity of the single particle level excitation energy. 
The reduced width for the single particle in any given nuclear state 
thus depends on the amount to which the single particle wave function 
is mixed into the state. The smaller the spacing between the 
nuclear levels, the larger the number of levels competing for the 
single particle wave function, and the smaller the reduced width 
to be expected for the single particle in any given level. If D 
is the energy spacing between particles of the same spin and parity, 
then y 2 /D is called the nuclear strength function for the interaction 
of the given individual particle with the nucleus. We expect the 
strength function to be large near the natural position of the single 
particle levels in the nucleus, and to be small far away from the 
natural position. If we do not know where the natural position of 
the single particle levels are, then we can assume that the strength 
function has been uniformly spread out over all excitation energies. 
In this case it has the value (Blatt and Weisskopf 1952) : 

y 2 2 x 10“^ 4 cm. (5.3.7) 



Of course one cannot expect either of these methods of 
estimating the strength function necessarily to give a good value 
of the reduced width for an individual level, since there are large 
differences in the degree to which the single particle wave function 
will mix into the wave functions for different states. Porter and 
Thomas (1956) have argued on reasonable grounds that the distribution 
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of reduced widths vi 11 be given approximately by the function: 



P(x) = exp - l/g x 
( 27 fx) 1 / 2 



(5 .5. a)* 



where x = 

This distribution function indicates that most reduced widths are 
considerably smaller than the average, and a few are much larger 
than the average. Such a great variability of individual reduced 
widths should always be borne in mind when estimating unknown physical 
parameters for astrophysical purposes. One should certainly not draw 
any strong conclusions which depend upon one or more reduced widths 
being close to the average value. It should also be emphasized that 
the above considerations apply to incident or outgoing particles which 
can be absorbed or emitted in only a single channel, which means that 
there is only one state into which the nucleus and particle can 
separate. When a particle can be emitted into many channels then 
the sum of the reduced widths for each channel starts to approach ' 
the sum of the average reduced widths. 



i/Ye now give for convenience the contribution that a distant 
level will make to the nonresonant parameter S, assuming -chat there 
is no interference between levels and that the contributions of all 
the distant levels aad linearly. We assume as before that the cross 
section d(a,b) for the reaction is given by the single-level breit- 
V/igner formula. Then 



S * 2.84 x 10 18 qp 2L ZtZ p (2J+1) 

(2s+l) (21+1) 



va 2 r b 

(e - EpV -4 1/4 r 2 



e v barns 

i 



If there Is serious interference, then the quantity in the square 
brackets would have to be replaced by a similar quantity obtained 
from the Breit-Wigner many-level formula. In cases where the 
reduced widths y " a nd Yv 8 are not known, then we can get an estimate 
of the contribution to 3°by estimating them on statistical grounds. 
The above expression for S should be evaluated at the Gamow peak 
energy E m . 



Two farther correction factors to the thermonuclear reaction 
rates will be mentioned briefly. They are small and of importance 
only for very accurate work. 



Salpeter (1952) has considered the errors which were made 
In the Gaussian approximation to the Integrand of (5,1.1), He finds 
that thermonuclear reaction rates should be multiplied by a correction 
factor: 
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P -7- = 1 + 5 - 35 + . . . 

12? 288V 2 



Corrections due to electron screening have been considered 
by ochatzman (1948, 1954a, 1954b), Keller (1953), and oalpeter (1954). 
Thermonuclear reaction rates are slightly increased because the 
penetration of Coulomb barriers is slightly assisted by the presence 
of electrons which neutralize the Coulomb field at large aistances 
from the nucleus. The formulation of the electron shielding 
corrections is complex and the resulting correction factors are 
usually not much greater than unity. No further consideration will 
be given here to these corrections. 

In many cases it is useful to be able to speak of the 
temperature sensitivity of a thermonuclear reaction. Let us write 

P = P 0 T n . 

Then the power n is obtained from the relation 

n = d log p . 
d log $ 

For nonresonant reactions it follows that: 

n = 1 { ?- 2). 

For resonant reactions the power takes the form: 

n = 0.0116 E r - 3 . 

' T 2 

6. HYDROGEN THERMONUCLEAR REACTIONS. 



< In this section we consider those hydrogen thermonuclear 
reactions which are of interest in the study of stellar interiors 
and of the changes in chemical abundances which can be produced by 
them. Those values of the nuclear reaction constants which appear 
to be the best available in the literature have been listed in 
Table 6,1. No attempt has been made to improve these constants 
except* in the case of the proton-proton reaction which is of 
particular importance in the construction of stellar models, and 
also in the case of the N 14 (p,y) 0^° reaction which is of similar 
importance. It turns out that none of the hydrogen thermonuclear 
reactions listed have any known nuclear resonances in the vicinity 
of the Gamow peak; hence the nonresonant reaction rate formula 
(5.1.2) applies to these reactions. This equation can be written 



TABLE 6.1 



Reaction 

i£(p,(3 + )d!* 

D^dlpjH* | 

Li 3 (p, a) He 3 
Li 7 (p,a)He 4 
Be 9 (p,d)Be 8 7 

B ?oi p,a ! LJ V 

BjJ(p,a)Be 8 

c }|jp»r)Ni3 

C i5! p » y ! N iJ 

N 1 ^ p, Y 015 - 

N} 8 p,a)C^ 

■J® P'T °}S 
oifip.rP^ 

0?-®(p, a )N 15 
Ne 20 (p, y)Na21 

s:Sfe?i;« , 

SiB&iiS# 

He^(a,v)Ll 7 
He 3 (He 3 , 2p) He 4 



S(e. 


v. barns) 




A 


B 


Note 


Reference 


3.12 


! X 10 -19 


3.3 


X 


10-13 


34.7 


1 


1 


7.8 


x 10“ 2 


7 


X 


10 4 


37.2 




2 


9.6 


x 10 4 


3.6 


X 


10 10 


42.6 


2 


2 


5 


x 10® 


6 


X 


!0l 2 


84.1 




3 


1.0 


x 10 5 


1.2 


X 


10 11 


84.7 




3 


3 


x 10 7 


3 


X 


1013 


103.6 




3 


2 


x 10 7 


2.5 


X 


10 13 


120.6 




3 


1 


x 10 8 


1.2 


X 


10 9 4 


121.0 




3 


1.2 


x 10 3 


1.7 


X 


10® 


136.5 




4 


6.1 


x 10 3 


9 


X 


10? 


136.9 




4 


3 


x 10 3 


5 


X 


109 


152.8 


3 


4 


1.1 


x 10® 


1.7 


X 


loj- 4 


153.1 




4 


1.1 


x 10 4 


1.7 


X 


1°1° 


153.1 




5 


y 


10 2 


>1.5 


X 


10 8 


166.7 


4 


3 


2 


x 10 5 


3 


X 


loll 


167 




3 


1 


x 10 8 


1.5 


X 


1014 


167 




3 


1 


x 10? 


1.5 


X 


10} 3 


181 




3 


1.2 


x 10 4 


1.8 


X 


1°1° 


194 


5 


6 


■5.6 


x lOf 


>1 


X 




194 




6 


5.5 


x 10 3 


^1 


X 


1°1^ 


194 




6 


3 


x 10® 


3 


X 


10 i2 


208 




6 


^10? 


•^1013 


208 




6 




0.6 . 


1.7 


X 


10® 


128.0 


6 


1 


1 


x 10 6 


4 


X 


loll 


122.5 


7 


1 
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Notes to Table 6.1 



1. The basic proton-proton reaction rate has been reevaluated as 
described in the text (Salpeter's result has been multiplied 
by a factor 0.72). The error in the reaction rate is now 
estimated to be 9 per cent, compounded from a 5 per cent error 




